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Stochastic effects during inflation can be addressed by averaging the quantum inflaton field over 
Hubble-patch sized domains. The averaged field then obeys a Langevin-type equation into which 
short-scale fluctuations enter as a noise term. We solve the Langevin equation for a inflaton field with 
Dirac Born Infeld (DBI) kinetic term perturbatively in the noise and use the result to determine 
the field value's Probability Density Function (PDF). In this calculation, both the shape of the 
potential and the warp factor are arbitrary functions, and the PDF is obtained with and without 
volume effects due to the finite size of the averaging domain. DBI kinetic terms typically arise in 
string-inspired inflationary scenarios in which the scalar field is associated with some distance within 
the (compact) extra dimensions. The inflaton's accessible range of fleld values therefore is limited 
because of the extra dimensions' flnite size. We argue that in a consistent stochastic approach 
the distance-inflaton's PDF must vanish for geometrically forbidden field values. We propose to 
implement these extra-dimensional spatial restrictions into the PDF by installing absorbing (or 
refiecting) walls at the respective boundaries in field space. As a toy model, we consider a DBI 
infiaton between two absorbing walls and use the method of images to determine its most general 
PDF. The resulting PDF is studied in detail for the example of a quartic warp factor and a chaotic 
infiaton potential. The presence of the walls is shown to affect the infiaton trajectory for a given 
set of parameters. 

PACS numbers; 98.80.Cq, 98.80.Jk, 98.80.Qc, 98.70.Vc 



I. INTRODUCTION 



Putting the successful inflationary scenario on the firm 
footing of a fundamental theory is one of the remaining 
challenges in cosmology. Recent years have seen consid- 
erable progress towards this goal with the construction 
of several concrete string inflation models; for recent re- 
views, see e.g. [H-Q and references therein. A top-level 
distinction among these models is the either closed or 
open string mode character of the inflaton field. A typi- 
cal example of the second class are brane inflation scenar- 
ios [^-01 '■ the inflaton field 4> corresponds (up to renormal- 
ization) to the distance between two branes embedded 
in a higher-dimensional background. While there is an 
ongoing debate about the form of the inflationary poten- 
tial v{4)) [Ml, a generic feature of these models is the 
field's kinetic term, which is of Dirac Born Infeld (DBI) 
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type rather than canonic [ij, |15| . One can understand 
these DBI dynamics as a geometry-imposed upper limit 
on the field's velocity 0, characterized by the so-called 
warp factor T((j)). This relativistic speed limit acts like a 
brake on the inflaton, forcing (f> to "slow-roll" even in re- 
gions where the potential is not flat. Hence, open string 
mode inflaton models provide an additional mechanism 
to generate quasi-exponential expansion. 

It is interesting to investigate the DBI analogues of 
standard inflationary calculations, such as the field per- 
turbations' evolution and spectra flpl [T7| . but also the 
effects of stochastic inflation |18l - |27| . Stochastic infla- 
tion provides a technique to assess quantum effects on 
the inflaton (j)'s trajectory averaged over a scale beyond 
the Hubble patch. In this way, one can define a coarse- 
grained field (p which consists exclusively of the large- 
scale Fourier components. To zeroth order, if obeys the 
slow-roll Klein Gordon equation, but its full evolution 
is subject to stochastic noise ^ from small-scale Fourier 
modes. For standard inflation with arbitrary potentials, 
this equation was solved perturbatively up to ©(C^) in 
Ref. [13]. Using this perturbative solution, Ref. [20l | 
showed how to obtain the Probability Density Function 
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(PDF) for ((/j) in the Gaussian approximation, with and 
without the volume effects {i.e. the size of each averag- 
ing domain) taken into account. The rehabihty of this 
treatment was further studied in Ref. [2]| . In this paper, 
we generahze both of these resuhs to the case of DBI 
inflation with arbitrary warp factors and potentials. 

The Langevin equation for DBI models was given for 
the first time in Ref. [28] , which aimed at studying "eter- 
nal inflation" p9l.[30] in the brane inflation context. Since 
then, stochastic effects in DBI (or general fc-inflation) 
models have been studied by several authors [311 ^32] . In 
this paper, we use the approach of Refs. [20l.l2l| to calcu- 
late the PDFs of a DBI inflaton field (with and without 
volume effects) for arbitrary functional form of the po- 
tential and warp factor. We illustrate our results by ap- 
plying them to the example model of "chaotic Klebanov 
Strassler (CKS) inflation", where T{(j)) and V{(f>) are 
known functions. The integrals encountered are, within 
certain limits, exactly calculable, and we discuss the be- 
havior and reliability of the resulting DBI inflaton prob- 
ability densities. 

However, we argue that all PDFs used so far in the 
literature suffer from a serious problem: they predict 
a non-vanishing probability for the moving I?3-brane to 
find itself outside the so-called "Klebanov Strassler (KS) 
throat" , i. e. outside the (part of the) extra-dimensional 
geometry with warp factor T{(j)). In fact, this problem 
is twofold. Firstly, at the bottom end of the throat this 
means that there is a non-vanishing probability to lit- 
erally find the brane "outside" the extra dimensions, in 
other words "out of space" , which is clearly meaning- 
less. Secondly, since the metric of the 6d bulk space [and 
hence the continuation of T{(j))\ is typically unknown be- 
yond the KS region, a string-inflationary scenario based 
on the brane's motion inside the throat becomes incon- 
sistent beyond its top end. Note that, while the latter 
question is of a more technical nature and we may hope 
to resolve it as our understanding of string geometries 
improves, the former issue is rather severe as the infla- 
ton's PDF does not respect the fundamentally geometric 
origin of the scenario. 

Put a different way, the compact character of the 
stringy extra dimensions (for the purposes of concrete 
model building, this means the well-known KS "corner" 
of the 6d manifold) can be translated directly into a re- 
stricted field range for the DBI inflaton. For consistency, 
any modifications of the classical trajectory induced by 
stochastic effects should still respect these geometry- 
imposed boundaries in field space. Hence, studies of 
stochastic DBI infiation so far were missing a tool to 
ensure the consequences of a stringy inflaton's geometric 
interpretation at the effective field theory level. 

To amend this problem, we propose to install "walls" 
at the boundaries {i.e. the bottom rg and the edge r^y) 
of the inflationary KS throat (whose radial coordinate 
is denoted by r). As a consequence, the stochastically 
corrected inflaton field value should remain within its 



the calculation of a new PDF respecting the boundary 
conditions imposed by the presence of the reflecting or 
absorbing walls. For the flrst time, we then determine 
this PDF in the presence of two absorbing walls at (/)o and 
using the method of images, see e.g. [H, HI]- As a 
by-product, we obtain the modified stochastic trajectory 
of the mobile brane within the KS throat and show that, 
in some cases, the presence of the walls has a significant 
effect. 

This paper is organized as follows. In the next Section, 
we start from the underlying background equations of 
DBI inflation and discuss how they can be used to formu- 
late their stochastic counterpart, i.e. the DBI Langevin 
equation. We pay special attention to the normaliza- 
tion factor of the noise term. We then solve the DBI 
Langevin equation up to second order in the noise and 
calculate the corresponding PDF along with its volume 
correction. In Sec. IIIIl we use these expressions to calcu- 
late the stochastic effects in the case of chaotic Klebanov 
Strassler infiation, both for a potential with and with- 
out a constant term. Where applicable, we calculate the 
domain of validity of the perturbative approach and con- 
sider the existence of a regime of eternal infiation [2^, [s^] ■ 
In Sec. IIVI we use the additional information on the in- 
fiaton's geometric bounds as an argument to implement 
two absorbing walls into the calculation of the PDF. In 
this way, it can be assured that even quantum effects do 
not violate these field space limits. Finally, we summa- 
rize our main findings in Sec. |Vl commenting on eternal 
inflation as well as on an overall picture for a realistic 
brane trajectory across the entire compact 6d geometry. 



II. DBI LANGEVIN EQUATION 

A. DBI Background Equations 

As a first step, we quickly recall the basic equations of 
DBI inflation at the classical level. Since models of this 
kind descend from {e.g. type IIB) string theory, they are 
originally represented by a higher-dimensional {d = 10) 
action for the stringy background and the embedded 
branes. After compactification to four dimensions, the 
effective field theory typically contains a gravity sector 
described by General Relativity and a four-dimensional 
inflaton field cf) {x,t) corresponding to the inter-brane dis- 
tance along one of the compactified dimensions. The 
model's effective 4d action therefore reads 



R 
2k 



v{cj)) ~ r(0) 



+T{^)Jl + 



T{<P) 



gt^'^d^ 



(1) 



where R is the four-dimensional scalar curvature and 



87r/m2 



being the Planck mass. It is useful to 



allowed range (j)o < (j) < (t>^^ at all times. This requires introduce (making use of the analogy with Special Rel- 
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ativity) the so-called "Lorentz factor" 7(0, d^fj)), defined factor 7 is easily expressed as a function of (j) only [28 



as 



vity 

m 



7(0, 9^,(1)) 



^1 + g^^-d^ </.a.0/r(</)) 



(2) 



Roughly speaking, this Lorentz factor measures how close 
the inflaton field's velocity is to the geometry-imposed 
speed limit ^JT(4>) (see below). In terms of 7 (the argu- 
ments of which we frequently suppress below), one can 
rewrite Eq. ([ij in the simple form 



7-1 
7 



n<t>) 



(3) 



Specifying to a Friedmann Lemaitre Robertson Walker 
(FLRW) universe with homogeneous scalar field matter, 
the Friedmann and Klein Gordon equations read 



7(0) = 



(9) 



Secondly, using the DBI slow-roll condition w 
kV{(I))/3, which is derived and justified in detail in Ap- 
pendixIXl one obtains the following first order differential 
equation for (p: 



r(0) 



(10) 



Obviously, except for the factor 7 appearing in the the 
denominator, this is the standard Klein Gordon equation 
in the slow-roll limit. A detailed derivation of Eq. ([TU)) is 
given in Appendix [Xj 



ry3 



-[(7-l)T(0) + m] , 



5H 



37-7' 



273 



■r'(0). 



(4) 
(5) 



where a prime denotes a derivative with respect to the 
inflaton field (p. Note that in the FLRW case, 7 from 
Eq. ^ simplifies to 
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1 - </)VT(^) 



(6) 



where the notion of y^T{(j)) as a speed limit is evident. 
Using the square root's expansion in Eq. (|6|), it is clear 
from Eq. ([5]) that, while 7 ss 1, the inflaton is close to 
standard dynamics. In this regime, obeys the usual 
Klein Gordon equation and hence can slow-roll only if 
the potential is sufficiently flat. Using the full expression 
(|6]), however, one can see that the inflaton velocity can 
never exceed y/Ticj)) even if the potential is steep. The 
limit where — >■ ^JT{(f) (and hence 7 — ^ 00) can there- 
fore be thought of as an additional, "ultra-relativistic" 
regime of inflation. Let us make this statement more pre- 
cise by taking the time derivative of the DBI Friedmann 
equation (|4]): the DBI condition to maintain accelerated 
expansion reads 



a 6 67 



(7) 



From this expression it is evident that the potential V{(t)) 
still has to dominate the energy density. 

Our next step is to notice that from combining Eqs. ^ 
and (O one obtains 



2H' 



(8) 



This formula has two important consequences. Firstly, 
replacing by this expression in Eq. the Lorentz 



B. Stochastic DBI Inflation 

We now proceed to applying the stochastic approach of 
Ref. [l^l to DBI inflation. In a first step, the "classical" 
inflaton fleld 0(i) is replaced by the coarse-grained held 
(p{t), which is a stochastic process. Based on the previ- 
ous considerations [see Eq. ([S])], we expect (p to obey a 
Langevin equation of the form 



2 H' 

K 7 



(11) 



where ^(i) is a noise term, describing the short wave- 
length part of the full inflaton fleld, and obeying the fol- 
lowing properties. 



(12) 



To proceed from here, the crucial question is how to de- 
termine the normalization factor C. 

In standard (non-DBI) inflation, this question may be 
resolved in several ways. One possible route to follow is 
to normalize the prefactor C according to the two-point 
correlation function of a massless test fleld in a de Sitter 
background. Indeed, in this case, the above Langevin 
equation (fTTjl (setting 7=1) reduces to <yj = C£_{t), with 
its trivial solution ip{t) — C ^(T)dr. Using Eqs. (|12p . 
one finds 



(13) 



Since the exact result is known and reads {(^^{t)^ = 
ff3t/(47r2) (setting ti„ = 0), we can read off immedi- 
ately that C ~ H'^/'^ / {2tt). Another way to see this is to 
interpret Eq. ([TT|) as a Brownian motion [2^, the field 
undergoing quantum kicks of amplitude H / {21:) in every 
Hubble time interval i?^^. This leads to the expression 
(^(p'^{t)'j ~ [i//(27r)]^ n, with the number of steps n given 
by n = t/ff-i = Ht. 



4 



In the case at hand here, i.e. stochastic inflation with 
a DBI inflaton field, we need to reproduce the same con- 
siderations, but for the modified Langevin equation (fTTj) 
that now comprises the Lorentz factor 7. For this pur- 



pose, let us remark that the fully covariant formulation of 
the DBI Klein Gordon equation derived from the action 
Eq. (m reads 



S^'^V^ - ^ g'^^g'^'^ (Vo V^</>) - ^ + |^ (7' -37 + 2) = 0. (14) 

To see that this gives back Eq. (O for spatially homogeneous DBI field (f) ~ (j>{t) in a FLRW universe, use the definition 
^ [or ([B]), respectively] of the Lorentz factor 7 and note that 

^'"^(VaV^^) v^^v,<^ = ^^r0. (15) 

At the perturbed level, setting 

(j){x, t) = (t){t) + 5(t>{x, t) , ds^ = - (1 + 2$) dt^ -I- a^{t) (1 - 2$) %da;Ma;^' , (16) 
one can show through repeated use of the Einstein equations for a DBI scalar field (see Appendix [B]) that 

J^fe + 3i/ (1 + h) 5^k + (^J^ + 2ei - ^ £2 - 2e? - I + ^ eiea - ^ + ^ + ^ <5iei - ^i^s 
+ + ^) = ^ (2 - 2ei + 62 + 2^0 (17) 

I 



in Fourier space with comoving wavenumber k. Here, we 
have used the {ei,Si) parameters defined in Appendix 1X1 
to write Eq. ([TT]) in a compact form. Note that in the 
limit where 7 — ^ 1, 5^ — 0, this gives back precisely the 
same equation as in the standard case (see Appendix |B|) . 

Ignoring the metric perturbation amounts to drop- 
ping the last term in Eq. (|17p . Moreover, in the "DBI 
slow-roll" regime, all the {€i,5i) parameters are small. 
Hence, if in addition we neglect all terms of at least lin- 
ear order in these parameters in Eq. (jl7p . we obtain the 
limit 

S4>k + SHS^k + ^Sq^k-0, (18) 

which again precisely corresponds to the standard equa- 
tion up to the replacement k fc/7. 

Usually, however, the perturbed Klein Gordon equa- 
tion is not written as in Eq. (1171) . but in terms of the 
Mukhanov-Sasaki variable Vk, which in the DBI case is 
defined as the combination 

Vk = aj^/^ S(t)k + z<^>k, (19) 

where z = aj^/ei. It can be shown (see Appendix |B|) by 
inserting this definition into Eg. (flTl) that Vk satisfies the 
following equation of motion |35| . 

dW ^ f^_}_^\ _o (20) 



Here 77 denotes the conformal time with which the scale 
factor is expressed as 0(77) — —l/{Hri) during exponen- 
tial infiation. In this regime we find d^z/drj^ = 2/77^, so 
that the normalized solution of Eq. ((20)) reads 

vUv) - -iL= (1 e-^'^-" , (21) 

where we set Cg = I/7 (which corresponds to the per- 
turbations' sound speed) and made the usual adiabatic 
choice of initial conditions. [Note that, as discussed in 
Ref. [l3], this choice corresponds to somewhat more re- 
strictive conditions for scalar DBI perturbations evolving 
according to Eq. ([20]) . as it would be the case for their 
standard counterparts.] Thus |5(/)fc|^ behaves as 

(22) 

in the long wavelength limit. Note that, unlike in the 
standard case, in the DBI picture the boundary between 
long and short wavelength regimes is not given by the 
Hubble radius but by the "sound horizon" CgH^^ = 
{jH)~^, as can be seen in the solution (PT|) . 

Now, the separation of the scalar field into long and 
short wavelength components. 



5 



(27r)3/2 



Q{k — sajH) ak5(j)k{t)e 



-ik-x 



ik-x 



(23) 



is precisely the essence of the stochastic inflation approach, where e in Eq. (l23l) is a small parameter (not to be 
confused with the first slow-roll parameter), and ak and are annihilation and creation operators, respectively. 
Using this expression we find that the noise term in Eq. pT|) can be expressed as 



(27r)3/2 



S{k — ea'jH) afc(50fc(i)e 



ik-: 



(24) 



with (50fe given by the solution of Vk, see Eq. (|2f p . Here we have used the condition that time dependence of Cs — l/j 
must be weak, |cs/cs| <C H (which corresponds to Si <C 1), which is required to justify our analysis of quantum 
fluctuations based on the mode function Eq. ((2T|) . We then find the correlation function of the stochastic noise is 
given by 



3 \SM' 



k—ea-yH 



(25) 



Thus we conclude that the analogue of the Langevin 
equation including a noise term in the DBI case is 



2 H' 

K 7 



27r 



(26) 



One notices that the factor 7 only appears in the classical 
term and not in the normalization of the noise term. 

By means of this equation, one can estimate in which 
regime the quantum effects are important. If the field's 
behavior is dominated by quantum effects, one can ne- 
glect the classical drift in Eq. (l26l) . As already men- 
tioned, using the properties ((T^ this leads to (^(p'^{t)) = 
H^t/{A'K^). However, the typical time scale /Si now is 
l/{Hj), since the "horizon" felt by the scalar field (and 
its perturbations) shrinks by a factor of 1 /j compared to 
the standard case. Therefore, the typical quantum kick 
undergone by the field in the characteristic time scale is 



A0q„ = 



m 1 



H 



47r2 jH 27r7i/2 



(27) 



On the other hand, if the quantum effects are negligible, 
then the equation determining the behavior of the field 
is nothing but the slow-roll equation of motion and this 
implies that 



V 



V 



(28) 



Setting A(/)qu — A(/)ci, one easily concludes that the cor- 
responding vacuum expectation value (vev) of the field 
4><, obeys the equation 



I 



4 V{(p.) [7(0,)]3/2 



(29) 



(Again, in the standard case 7=1 and one recovers the 
usual criterion.) Eq. ([^^ allows us to decide for which 



values of the infiaton quantum effects play an important 
role. 

So far we have used the cosmic time t as the time vari- 
able in the Langevin equation, see Eq. (pS)) . The choice 
of the time variable, however, is a subtle issue as differ- 
ent choices may lead to physically inequivalent results. 
It has recently been advocated in Refs. [l^, [l^l that in 
many cases it would be more appropriate to use the num- 
ber of e-folds, = Ina, as the time variable depending 
on what quantities we wish to calculate. Indeed, in this 
case, the results obtained from the stochastic formalism 
coincide with those derived from perturbative quantum 
field theory [2^ [23| . Written in terms of the number of 
e-folds, the Langevin equation ([26)) reads 



d(/3 
diV 



2 H' 
1^~H 



H 
2^ 



(30) 



where ^(A^) is a new stochastic process (for which, al- 
lowing for slightly slippery notation, we still denote by 
the same symbol) such that {iiN)^{N')) = S{N - A'). 
It is easy to check that a free field satisfies ((p2(7v)) = 
ff2(iV - Ain)/(47r2) as expected. Below, we will carry 
out our calculations in both time variables, cosmic time 
t and the number of e-folds N. 



C. Solving the DBI Langevin Equation 

We now solve Eq. using a perturbative expansion 
in the noise as shown in Ref. j2^]. We use the following 
ansatz for the Hubble-patch averaged field. 



(fit) = (pcl(t) + Sipi{t) + S<fi2{t) + 



(31) 



where (ficiit) is the classical field, Sipi{t) cx 0{£,) and 
5(p2{t) cx 0{£,^). In principle, this expansion can be car- 
ried to any order in ^. At zeroth order, we get back 
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the classical slow-roll equation ([T0| . At first order, one 
obtains an equation for dipi(t), namely 



dSipiit) , 2i7'(<^ci) 



dt 



K 7(Vci) 
27r 



i?'(<^cl) 



7(</'c 



(32) 



As expected, this equation differs from the standard one 
by the presence of the Lorentz factor and its derivative. 
At second order in the noise, one obtains the equation 
describing the evolution of 5ip2{t), 



dt 



K 7(Vcl 

3 

An 

1 H'{^,i) 



iJ"(^el) 7'(<^cl) 



H'i(Pc\) 7(</'ci) 



K 7((/?cl) 



,7''(^ci) 



H'((pci) lifci) 7^(<^ci) 



(33) 



The same remark as before is valid: the equation for 
6(p2{t) is modified by the presence of the factor 7 and its 
derivatives. Clearly, the equation contains derivatives of 
7 up to second order because it is second order in the 
noise expansion. 



We are now in a position to solve the above equations. 
Since they are first order differential equations, they can 
be solved by varying the integration constant. One finds 
for 6ipi{t) that 



dt ^TTT 7['^cl(i )] C(t ) 



(34) 



while for 6(p2{t) one obtains 
3 H'[^,i{t)] 



4vr "f[iPciit)] Jt,^ 



dt'H'/^[v,iit')]j[ip,,{t')]at')SMt') 



l^l[fc\{t)] 



dt' 



H"'[v,i{t')] 7"bci(t')] 



H'iipMt')] liMt')] 7'bci(t')] H'iMt')] 7[Vci(i')] 



S^iit'). (35) 



r 



From the above expressions and the properties of the 
noise given by Eq. (IT^ . it is obvious that {6lpi) — 0, and 
for the second moment we find 



(36) 



The next step is to calculate {6(p2)- Lengthy but straight- 



forward calculations lead to 
H' 



{S1P2) = 



7 



dV- 

{H'/l)' 
2{H'h) 

H'h 
Anm? 



H' 




7 


) /: 


f( 
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3 gV 
2 H' 



{5^1) 



H'2 



(37) 



As in the standard case, everything can be reduced to 
the calculation of a single quadrature. As expected in 
the DBI case, this quadrature contains the factor 7. 

Using these results, one can now calculate the PDF, 
Pc{(p, t), which describes the probability of the stochastic 



7 



process to take a given value ip at a given time t in 
a single coarse-grained domain (see Ref. |2g|). 



1 



exp 



{5^2)). 



(38) 



[In the last line, we have introduced the definition Pg 
for later use, see Eq. ((5^ .] If, however, one is inter- 
ested in spatial averaging over the entire Universe (in- 
stead of a single domain), one has to include a weight 
factor o?{ip) = exp [3 / dr iJ((p)] for the physical volume 
of each Hubble-sized domain. This leads to 



J(^-(^[e])e3/d-^^(vKl)\ 



^g3/dr//(¥.K])^ 



: exp 



(^_(^)_3JTJ) 



2n 



v/27r {5^1 
Pg(^-((^)-3/Tj)^ 



2(M) 



(39) 



where again the last line is a definition used in later Sec- 
tions. In Eq. {if) is the single domain averaged 
mean value (i^) = (^ci + {^'^2)1 to which the volume ef- 
fects induce the additional correction given by 



2,1^ J 



3 j dTH'{T){5^^{t)5ip^{T)) 



127r7iJ 



-(5v'?)-(40) 



Once this integration carried out, the volume-weighted 
distribution Pv{(p,t) as defined in Eq. (p9)) is also com- 
pletely determined. Again, comparing the results of this 
section with those for the standard case presented in 
Refs. [23,1m, see that the only changes are the addi- 
tional powers of 7 found in Sipi(t) and Sip2{t). Below, we 
calculate Pc{(p,t) and Pv{(p,t) for an exemplary shape of 
V{(p) and T{ip). 

As discussed above, the Langevin equation can also be 
written with the number of e-folds as the time variable. 
It is straightforward to repeat the above analysis for the 
corresponding Langevin equation pop . In particular, the 
first and second order corrections obtained in terms of 
e-folds read 



SifiiN) 



1 H' 



N 



AN 



,H \N')^{N') 
H'{N') 



m'). (41) 



and 



1 H' 



dN' Hj6ipi{N')£_{N') (42) 



dN 

A'in 



In Eq. (l42l) . arguments in the integrands have been par- 
tially suppressed where they are evident. As before, it is 
easy to calculate {6ipl) and {S(p2) from these expressions. 
It is found that 



N 



87r2 i72^2 



H'3 



dV'. 



(43) 



For the case of a standard kinetic term (hence, 7 = 1), 
this equation coincides with Eq. (40) of Ref. ^7^. Fur- 
thermore, we find from Eq. (|^^ that 



N 



H' 



327r2 



i/'2 



(44) 



Again, if 7 = 1, this equation corresponds to Eq. (48) of 
Ref. [23. Let us notice that, as in our previous results 
of Eqs. ([36]) and ([37| . the expression p3t for (^5lp\) ^ is 
sufficient to obtain {Sip2)j^j in Eq. (j44p . i.e. no additional 
quadrature is necessary. 

Following the same steps, one can also evaluate the 
PDFs in terms of e-folds. The corresponding probabil- 
ity Pc^\ip,t) is similar to Eq. (|38l) . except that now the 
formulas (|43|) and ((44|) should be used in their respective 
places. The definition of the volume weighted distribu- 
tion also remains the same, see Eq. (l39l) . but the term 
3/""" J now reads 



(3/-J) 



AT 



12 H' /■'^'" , , 
— -. / dib — In 



127r 777 



ml. H 



- {6ipl) In 



fP3 

H 

771„, 



H 



(45) 



This equation should be compared to Eq. PO)) . Again, 
assessment of the volume effects requires the calculation 
of a new quadrature. 



III. APPLICATION TO BRANE INFLATION 

In the following, we apply the formalism developed 
in Sec. |TT] to a popular class of string-inspired inflation 
models with DBI kinetic term. We focus on scenarios 
of the brane inflation type, for which the inflaton field 
corresponds to the position of a £'3-brane embedded in a 
higher-dimensional background. Successful model build- 
ing requires that the six extra dimensions be deformed in 
a way described by a warp factor r(0); the resulting ge- 
ometry is commonly called a Klebanov Strassler throat. 
Scenarios of this type have been the subject of a vast 
body of literature, see e.g. Refs. [sl-fisl [28l. |36| - [39| . For 
our purposes, we denote the warp factor and the infla- 
tionary potential by 



A ' 



F(0) = 



m^(j)\ (46) 
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where e = ±1. The plus sign identifies so-caUed "Ultra- 
Violet" (UV) models (where the D3 moves from the edge 
towards the bottom of the throat geometry and hence the 
field value decreases during infiation), while the minus 
sign refers to the "Infra-Red" (IR) setting (where the 
D3 climbs out of the throat and the inflaton's field value 
grows with time). 

For completeness, the formulation Eq. P5)) of V{4)) 
includes a Coulomb potential term due to the I?3's at- 
traction towards a Z)3-brane sitting at the bottom of the 
throat. Originally, the small Coulombic attraction result- 
ing from this very fiat (oc 1/'/'^) potential was considered 
the infiaton's only driving force, ignoring the (potentially 
much steeper) second term cx 0^ in Eq. (j46| [3]. In fact, 
it was shown that in this case the DBI dynamics do not 
affect the infiationary evolution and the stochastic 
effects were also assessed in the same reference. There- 
fore, in the following, we ignore the oc 1/0^ Coulombic 
contribution in Eq. pS)) . Note, however, that conceptu- 
ally the presence of an anti-brane at the bottom of the 
KS throat shall be important for our reasoning in later 
Sections. The quadratic potential term in Eq. (^5)) has a 
different status. We shall treat it here as a phenomeno- 
logical description of the potential that the various back- 
ground moduli fields produce for the mobile D3-brane. 
The exact shape of these moduli contributions is still a 
subject of active research |7hl3i] . 

From the cosmological point of view, the infiation 
model of Eq. (j46|) has three parameters, the mass m, the 
dimensionless constant A and Vq (with dimension WpJ- 
In fact, as we show below, it is rather two dimensionless 
combinations of (m. A, Vo)j which we shall call (a, /?), that 
characterize the evolution. We define these parameters 
by 



Unml, 967r2 



Am^ 



(47) 



Physically, /? measures the importance of the constant 
term relatively to the mass term in the potential (recall 
that we are neglecting the Coulomb term which involves 
the parameter /i). 

The geometric interpretation (in terms of an extra- 
dimensional brane position) of the infiaton field enforces 
some intrinsic consistency conditions. The renormaliza- 
tion relating the infiaton field to the (radial) throat 
coordinate r reads 4) = ^/T^r, where T3 = [(27r)'^(7sQ;'^]~^ 
is the tension of the D3-brane, calculated from the string 
coupling and scale a' . (We do not consider a possi- 
ble motion of the brane along angular coordinates of the 
throat.) Let the edge of the throat correspond to some 
r^jy , where the KS corner of the geometry is connected to 
the compactified six-dimensional bulk. Since the metric 
outside the throat is unknown, one has to impose r < r^^ 
to ensure that the brane stays inside the well-defined KS 
region. In terms of stringy background parameters, one 
can express 



where is a positive integer representing the total 
Ramond-Ramond (RR) charge and v represents the (di- 
mensionless) parameter measuring the volume of the five- 
dimensional submanifold that forms the basis of the 6d 
throat in units of the five-sphere volume. Via the infla- 
ton's renormalization, this evidently is an upper bound 
on the infiaton (p. Note that, depending on whether we 
are talking about UV or IR models, this constraint affects 



the initial or final field value, 



i(uv) 



< 



.(m) 
'^cnd 



< 



A second consistency condition is the requirement that 
the volume of the throat to be smaller than the total vol- 
ume of the compactified extra dimensions. This total 
volume has observational significance since it enters into 
the four-dimensional Planck mass mpj. From this con- 
straint it follows that 



< 



< 



(49) 



This means that inflation always occurs for sub- 
Planckian values of (p. On the other hand, the bot- 
tom of the throat being located at rg, one must have 
> 00 = 'Jr^fQ. Moreover, for the model to be 
valid, the (physical) distance between the brane must 
be larger than the string length and one can show that 
this amounts to 



> 



(50) 



Note also that the parameters of the warp factor and 
potential in Eqs. (j46|) can be calculated in terms of the 
stringy parameters. Physically, T{(j)) is the position- 
dependent brane tension and it can be written as T((/)) = 
T3 {<j)/(j)ijv) , which imphes [compare Eqs. p6)) and 
that 



A = 



(51) 



The constant term Vb is given by Vo = Att^v^q/N , which 
can also be expressed as 



yo = 2;i4(ro)T3 



(52) 



where hijp) = (f)/(j)^^ is the warping function as it appears 
in the lOd metric [it holds that T(0) oc h'^icj))]. 

It turns out that Eq. (|49|) can be rewritten in a 
more quantitative way. Since the volume of a Klebanov 
Strassler throat is known, yg*hroat _ 2Tr'^g^Afa'^r'^^, and 
the total six-dimensional volume is related to the four- 
dimensional Planck mass, Vq°^ = m,'^^{2ny g'^a'^ / {IGtt) , 
one deduces for our parameters (a, /3) deflned in Eq. (j47p 
that 



1 tf(.-o) 



a V247r3 Af 



< 1 



(53) 



(48) 

A similar equation [see Eq. (2.10)] was used in Ref. [3£ 



9 



Chaotic Klebanov Strassler Inflation 



and the final result for (Sipf) is 



In our first example, we set the parameter 13 — 0, hence 
we are considering the case of "pure" CKS inflation with- 
out a constant term. The potential and warp factor hence 
are given by 



(54) 



In the slow-roll limit, it follows from Eq. (|4]) that ~ 
ATTm'^cj)'^ / {^rn^J (see Appendix [X| and, therefore, the 
Lorentz factor calculated from Eq. ^ behaves as 



2Xm? 1 



3k 



(55) 



If we plug these expressions into Eq. (p6|) . the calculation 
of (Sipi) can be reduced to an exactly solvable integral. 



4m 
3 



2 r 



-In'^' 

2a \ ifci 



Aa V 1 + 7 



(56) 



where 7 = j{ipc\) and = For {^Lp-i) we have 

from Eq. ^ that 



7^-1 1 
7^ Vc\ 



37TOp 



(7fn'^f„ - 7Vcl) 



To calculate the volume effects, another integration is 
necessary. With the potential and warp factor given by 
Eq. (ISH), one finds from Eq. (gOl) that 



31^ J 



7 m? 
-B 





1 mpi 

\ 'Pin 


A. 


3 3\ 


PI 


'4' 2/ 



3/2 



TOpj 



-127r7^-^^^ei, 



3/2 



- ^(-l)-V4«-5/4 



3 3 
"'4' 2 



(58) 



where B is the incomplete Euler's integral of the first 
kind defined by B(z,a,6) = ^"-^{1 - tf-^dt HI Si- 
Notice that, in the above equation, the function B takes 
in fact complex values. However, multiplied by the factor 
(— 1)"^/'^, the result is real as it should be. 

Let us discuss these results in more detail in the light 
of the consistency constraint Eq. (|49|) . In the limit of 
small field values compared to the Planck mass, the 
Lorentz factor is large and can be approximated by 
7(cpci) — "T-pi/ (v^'Pci)- It hence easy to show that 



\mp, 



3m2 



(59) 
(60) 



The expression giving the volume effects can also be sim- 
plified. Notice that B{z, a, b) — z°'2Fi{a, 1 — h^a + 1, z), 
where 2-F1 is the hypergeometric function [4ll,|4l. Using 
the asymptotic behavior of 2F1 [lllHl, straightforward 
calculations show that 



a 



(61) 



Let us now establish the results in terms of the numbers 
of e-folds. Using Eq. one obtains 



N 



2 


m"^ TOpj 


U3 


^^1 


7?„ ^fn 


7 fli 


3 


7^ (^2 


3 


PI 


3 m6 

PI 
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a 

PI 



1 



v3/2 



arcsinh 



ai/2 



a 



3/2 



arcsinh 



1 TO 



a 



1/2 



(62) 



One can easily check that this expression vanishes at ini- 
tial time as expected. The corresponding expression for 
{8(^2) can be deduced from Eq. (pi)) . The result reads 



Vn 



1 7^ 



(Spl) 



4- (7.14 -7^4) (63) 



In the limit where the Lorentz factor 7 is large, the above 
expressions can be approximated by 



2m^ 



{S<f2) 



N 



2w?- ipc\ 
3toJ, Ja 



In 



Pc\ 



ln(^ 



(64) 
(65) 
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It is interesting to compare these formulas to Eqs. ((59)) 
and ([50)1 . We see that working in terms of the num- 
ber of e-folds simply introduces (apart from numerical 
prefactors) a logarithmic correction to the correlation 
functions: roughly speaking, the new correlation func- 
tions are obtained from the old ones with the replacement 
LPci — > (/^ci ln('^ci/</5in)- This is confirmed by a calculation 
of the volume effect. Using Eq. (|45p . one obtains 



167r 

2" "^t:! 



X <! In ( ^ I In 



A \ 1/2 

47r \ m(/Jci 



2 



47r \ miyJi, 



^ \ 1/2 



(66) 



This formula should be compared to Eq. (|5T|) . 

In order to see whether the stochastic effects are impor- 
tant or not, we must normalize the model's parameters 
to the COBE observations. This was done in Ref. [l7| . 
where it was shown that [see that reference's Eq. (127)] 



TO 
TO 



PI 



45 1 
4^r2 ~ 

CMB 



(67) 



where the quantity Q"^ /T'^^-^ can be expressed in terms 
of the CMB quadrupole according to 



Q 



T 



^^6x10-^ 
47r 



(68) 



In Ref. [itS [see Eq. (116)], it was also demonstrated that 
the first slow-roll parameter ti (see Appendix Elfor a pre- 
cise definition of the slow- roll hierarchy), for the model 
under consideration, can be expressed as ei ~ v^/(47r). 
Therefore, in a realistic inflationary situation we always 
have a <C 1. 

As a rule of thumb, the inflaton field can be said to 
behave quantum- mechanically if the correction ((5(^2) to 
the mean value is of the same order as this classical field 
(^ci- For a more detailed argument, let us consider the 
calculation carried out in Ref. [i^ (see also the discussion 
in Ref. [13]). There, the authors compute the number of 
e-folds {N) = J dt{H), which we adapt to the case of 
DBI inflation as 



, , {H} 
77r7ci • 



(69) 



In the present case, i.e. for the potential and warp factor 
of Eqs. this gives 



dV' 



m tpci 7ci 
H' 



1 



(70) 



If (^1^92) ^ 'Pel, then (iV) = Nc\ and the trajectory is 
indeed classical, confirming our rule of thumb. If we work 
out the above condition ignoring unimportant numerical 
factors, we find that 



£1 



7^2 

CMB 



< 1 



(71) 



We conclude that the stochastic effects do not play any 
important role in the model with /3 = and, therefore, 
that eternal inflation is not possible. This conclusion, al- 
though obtained with a different method, is in agreement 
with Ref. [i^]. In view of this result, we do not discuss 
the domain of validity of the perturbative approach in 
the (3 = case, but instead turn straight to the case of a 
CKS potential with constant term. 



B. Chaotic Klebanov Strassler Inflation with a 
Constant Term 

As a second example, we keep T{(j)) as in Eq. (j54p . but 
take /3 ^ 0, so that the potential is given by 



(72) 



With this potential, the Lorentz factor [compare Eq. ([55 
is given by 




1 



(0/TOp 



a 2/3 e (^/TOpJ-^ 



(73) 



Then, in order to determine the corrections to the vari- 
ance {Sfi) and mean value {Sip2), one has to compute 
the kernel of the integral in Eq. l\'SQ\i . One obtains 



ff7 

IF 




(74) 



Unfortunately, the (third power of this) expression is 
too complicated to perform the integral Eq. ([36)) ex- 
actly. However, while the constant term of the potential 
Eq. (|7^ dominates (/? ^ 1), it is legitimate to approxi- 
mate Vq + em'^(j>'^/2 by Vq. In this case, the calculation 
of Eq. (|36p can be done and leads to 



15 ^-07c^, ' 

1 + 2af3 {(pc\/m 
2aP {(fici/m^^] 



1 + 2af3 (yin/TOp,) 
2a/3 {(Pin/nip^ 



5/2 



(75) 
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In the same limit, one can also estimate {8^2) 
reads 



The result 



{5^2 



m 

TO 



PI 



TOpj 



fc]_ 

7ci 

1 

2^ 



TOp, 



7i 



(76) 



where the explicit expression of (S'-Pi) is given in Eq. (1751) . 
The evolution of {5'-p\) and {5ip2) in the IR and UV cases 
is represented in Figs. [1] and [H 

Finally, we have to estimate the volume effects. For 
this purpose, we have to compute the correction to the 
mean value of the field given by Eq. (j40[) . In our approx- 
imation, one has H^-f^/{H'f ~ y/nVo/?, [H^-i^/{H'fV 
This means that one can express the integral of Eq. (j40p 
in terms of the integral appearing in Eq. p6p. i.e. in 
terms of (^Sipf). After carrying out this calculation, one 
obtains 31^ J ~ 0. Therefore, at first order in our ap- 
proximation, volume effects simply are absent. In fact, 
one could have guessed this result from the very begin- 
ning: at first order, the Hubble parameter is given by 
H ~ ^kVq/S which is a constant. Looking at Eq. (l39l) . 
we see that, in this case, the weight related to the vol- 
ume can be taken outside the integrals in the numerator 
and denominator. Hence, the corresponding contribu- 
tions cancel out and, at this order, there is no volume 
effect. 

For the same reason, there is no difference between the 
correlation functions computed with the time variable t 
or with the e-fold variable N, see for instance Eqs. (|36|) 
and (|43|) . This illustrates the fact that the difference 
between the two approaches can be important only if the 
Hubble parameter evolves significantly. This is the case 
for chaotic inflation 0,[23] (although for the DBI version 
of the chaotic model studied in the last subsection, the 
corrections were only logarithmic) but not for the model 
under scrutiny here. 

Before investigating in more detail the previous result, 
let us establish its domain of validity. It was shown in 
Ref. I^T] that one can trust the perturbative stochastic 
treatment as long as the mean value of the infiaton field is 

such that {(fi) e [(fcl ~ |A(/?min iVd)] , + At^max (<Pcl)] 

where Aip^^i^ and Aip„ 
tions. The first of those reads 



^min Vr^ciyl ? V^cl r ^Y^max ' 

can be found from two condi- 



max 

a:6[¥'cl,¥'cl+Aip(¥'cl)] 



6 



TT/ir 



A^^(77) 



where H^^^^ denotes the fourth order derivative. The sec- 
ond condition can be expressed as 

1 " 



max 

a:6[¥'cl,Vcl+Aip(ipcl)] 



(78) 



These two conditions must be simultaneously satisfied 
and, therefore, the tightest bounds on Ai^Jmin and A(^max 
that follow from Eq. (fTT]) and Eq. (fTSt , respectively, give 
the reliability of the perturbative treatment. Replacing 
again Vb + £to^^^/2 by Vq where applicable, straightfor- 
ward manipulations show that Eq. ([77]) leads to A(^niax = 
A(^min = 3(^ci while Eq. ([78]) gives At^^ax = Aip^in = 
2ipci. Therefore, one concludes that the perturbative ap- 
proach is correct as long as (ip) e [tpd — 2(pci, fd + 2(/3ci]- 
The allowed region is represented by the hatched area in 
Figs. ID and H 

We now return to Eqs. ((75|) and ((76|) . These expres- 
sions can be further simplified if we take into account the 
fact that the vev of the inflaton field (measured in units 
of the Planck mass) must be small, compare Eq. (|49|. In 
this case, the Lorentz factor from Eq. ([75)) is given by 
7 ~ TOp|/(V2a/3(pci) > 1 and 



TO 
TO 



PI 



16 
I5V2 

16 / TO 

15V^ V'^p 



a 



1/2 



(79) 
(80) 



In the same way, one can also calculate the correction to 
the vev of the inflaton fleld. One obtains 



- ( 


\ 2 

TO \ 


^3/2 ^3^ ^3^ 


15\/2 V 


TOp J 
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( — ) 


2 ^3/2 ^3^ 


I5V2 




aV2 ^2^ • 



(81) 
(82) 



Notice in particular that, in the UV case, the correction 
is negative. This is conflrmed in Fig. [2] 

With the help of these approximations, let us estimate 
when the stochastic effects are important. As discussed 
before, one expects the quantum effects to play a role 
when {S1P2} /fci — 1- In order to use Eqs. (|8T|) and ((82)) 
for this purpose, we must again calculate the COBE nor- 
malization; this will allow us to re-write the parameter 
combination m^P^J^Ji 
and dM]). In Ref. [13 , 



^2^3/2^^^2^0,1/2-^ appearing in Eqs. ((5T)) 
In Ref. lia, it was shown that [see Eq. (153)] 



TOp 



45 

167r T2 

CMB 



/TOp^\2 a_ 



(83) 



where (p^, is the value of the inflaton field when scales 
of astrophysical interest crossed out of the DBI sound 
horizon. Moreover, it was also demonstrated that — 
1 ^ A61 [see Eq. (155)]. Therefore, we finally arrive at 



TO 
TO 



PI 



/33/2 457r3 g2 



a 



1/2 



2^2 

CMB 



V a 



3/2 



(84) 



As a consequence, for the IR case described by Eq. ([81 
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FIG. 1: Evolution of the (quantum) scalar field in the IR case (hence inflation proceeds from left to right) for a — 38, /? — 3.7 
and m ~ 2.19 X 10~^mp|, as implied by the COBE normalization. Note that the choice of the parameters is such that the 
condition /3/a <§; 1 is valid. For this example, the spectral index and the running of the model are respectively given by 
rig — 1 ~ —0.11 and Qg ~ 0.0023, as it can be shown using the results of Ref. [TJ]. The initial condition is ipin = 10~'*mp|. 
The green dotted line represents the classical evolution without the quantum effects. The red dashed line represents the mean 
value of the quantum scalar field, namely (pci + (fiip2) while the two blue dashed dotted lines on both side of the mean are 
ipci + {5^2} ± \/ (^Vi). On the right panel, the hatched region represents the region where the perturbative treatment used in 
this article is valid. For the parameters chosen here, the perturbative approach breaks down at <^ci — O.OOlmpj. 



the condition {Sip2)^^ /^Pci ^ 1 is equivalent to 

x(^) . (85) 

Let us discuss this expression in detail. It is interesting to 
note the dependence on the initial value of the field (p-m- 
The smaller (pi„ is, the sooner the influence of quantum 
effects sets in. (Recall that ipi^^ is bounded from below 
by the bottom of the throat, (p-m > 4>o, and that the 
inflaton field value in this scenario increases as the Di- 
brane climbs out of the throat.) For (y9i„ = 10~^'^mpj, 
corresponding to the left panel in Fig. (TJ one obtains 
(y^ci/Tipi ^ 0.008 in good agreement with the plot. For 
an exemplary initial field value of t/Jin = 10~'*mp, , one 
has ipci/rrip^ > 0.0005. Therefore, if the brane starts 
its evolution deep inside the throat, then the stochastic 
effects are dominant. 

For the UV case of Eq. ([5^ , one can repeat the same 
discussion. One finds that this case the limit is given by 

mpj Vv2/ VT'cmb/ V"/ 

(86) 

Contrary to the IR case, we see that there is no depen- 
dence on the initial conditions anymore. For the param- 
eters in Fig. [21 one obtains (pd ^ 2.19 x 10~^mp, in good 
agreement with the plot. In addition, one can also check, 
by comparing the two panels of Fig. [5J that this value 



does not change much when ipi^ is modified. The con- 
clusion is that in the UV case, the stochastic effects only 
play an important role when the brane is approaching 
the bottom of the throat, i.e. towards the end of brane 
inflation in its UV incarnation. 



IV. STOCHASTIC DBI INFLATION AND THE 
FINITE SIZE OF EXTRA DIMENSIONS 

We have established that, in the model of Eq. (|72l) . 
the stochastic effects can be dominant in a relevant and 
realistic inflationary regime. Our next goal is to compute 
how the (classical) behavior of the inflaton field is modi- 
fied in the presence of stochastic noise. In principle, the 
above considerations already answer this question since 
we have computed in detail the time evolution of {5ip2) 
and (Sipi). However, it is clear from Figs. [1] and [2] that 
a first issue is the limited validity of the perturbative 
regime. Outside the hatched region in the right-hand 
panels of Figs. [T] and [H one can no longer follow the 
evolution of the quantum field. In particular, it seems 
obvious that the regime of eternal infiation cannot be 
described in this framework, even if we have established 
earlier that it is likely to exist. However, even letting 
aside the issue of validity, there is another, much worse 
problem that renders our previous treatment highly un- 
satisfactory. The probability density functions given by 
Eqs. (1551) and can a priori extend into the range 
where Lp < (possibly even into the range (/? < 0!). 
Clearly, since the inflationary scenario at hand is built on 
the notion of the inflaton as the KS throat's renormalized 
radial coordinate rp < r < r-^^, a fleld value ip < (pQ is 
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FIG. 2: Evolution of the (quantum) scalar field in the UV case (hence inflation proceeds from right to l eft) for a = 38, /3 = 3.7 
and m ~ 2.19 x lO^^rripj as implied by the COBE normalization (this case corresponds to Fig. 9 of Ref. [IJ)- For this example, 
the spectral index and the running of the model are respectively given by rig — 1 ~ 0.11 and Og ~ —0.0023. The initial condition 
is (^in = 10~^mpj. The conventions are analogous to those in Fig. [T] The fact that {5if2)_^^ < 0, see Eq. (|82p . and, therefore 
(<^) < ifici, can be clearly observed in the above plot (red dashed line). On the right panel, the hatched region indicates the 
domain of validity of the perturbative approach. 



inconsistent. To our knowledge, all works on stochastic 
inflation in brane inflation published so far are subject 
to this issue. In fact, we face here a deep conceptional 
challenge: in brane inflation, the finite size of the extra 
dimensions plays a fundamental role. (In Ref. 40], for 
example, it was shown how the finite size condition can 
be exploited to cut down the allowed parameter space 
at the effective field theory level.) How can this crucial 
importance of the geometric restrictions be implemented 
into the description of stochastic DBI infiation? 

We propose to address this issue by introducing a (re- 
flecting or absorbing) wall at 0o, fohowing the reason- 
ing of Ref. [13]. (A similar technique has been used in 
Ref. [15 in order to study the quantum behavior of the 
quintessence field. Analogously, the approach we present 
below may be applied to the case of "DBI-essence" in- 
troduced in Ref. (4^1 ■) This wall at the bottom of the 
throat then marks the "end of the world" as imposed by 
the finite string geometry and prevents the stochastically 
corrected value of from reaching values smaller than 0o • 
Let us now explore the consequences of this proposal. 

Let us recall from Refs. (i^ . |45| that, if we start 
with a normalized distribution P{x) for the variable x, 
i.e. P{x)dx = 1, then the normalized probability 
density function in the presence of a refiecting wall at 
X = a is given by P{x) + P{2a — x), provided that 
P{2a — x) is also a solution of the relevant (approximate) 
Fokker Planck equation. (This is the case for a Gaus- 
sian PDF.) One can check explicitly that this distribution 

is correctly normalized, °° [P{x) + P{2a — x)]dx = 1. 
As a thought experiment, let us install a reflecting wall 
at the bottom of the KS throat at the position 0o. Then, 
according to Refs. [3 EEIj our new distribution would 



be given by 



-Pwall(V') = 



VMM) 



cxp 



■ exp 



(200 - - </?cl 



2 (5^2) 



2 {b^D 



.(87) 



Notice that the wall distribution is the same with or with- 
out the volume effects since we have shown before that 
3/^J~0. 

However, in the present context, when the brane 
reaches the bottom of the throat, it annihilates with 
an anti-brane fixed at 0o. Note that, while the small 
Coulombic contribution of the inter-brane potential in 
Eq. ()46|) was ignored in our calculations above, concep- 
tually the presence of the anti-Z33 at the bottom of the 
throat is important for our present argument: therefore, 
a correct physical description of this situation involves 
an absorbing wall rather than a refiecting one. A more 
difficult question concerns the boundary condition to be 
chosen at (^^^ (where the KS throat joins the 6d bulk 
with unknown metric). 

A more appropriate description of the situation could 
possibly be given by a "multi-throat" scenario: for in- 
stance, one can image that the mobile D3 brane trav- 
els through a 6d compact space made of a (negligibly) 
small unknown bulk and two well-defined KS throats of 
different depth. To account for the anti-£'3s present at 
the bottom of each throat, absorbing walls are installed 
at the respective positions (or their field value counter- 
parts, respectively). Clearly, the detailed treatment of 
such a situation will involve a high degree of technical 
sophistication because a continuous description of warp 
factor and inflaton potential is needed across the entire 
6d manifold. 
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FIG. 3: Evolution of the (quantum) scalar field in the UV case (hence inflation proceeds from right to left) for a = 38, /3 = 3.7 
and m ~ 2.19 x lO'^rripj as implied by the COBE normalization (this case corresponds to Fig. 9 of Ref. [IJ)- The conventions 
are analogous to those in Fig. [1] The black solid line represents (y) given by Eq. (|10ip in the case where the two walls are 
located at (jio = 10~^mpj and (j)^^ — 10~^. The initial condition is (jj-m = 5 x lO^'^mpj. This figure should be compared to 
Figs. [21 The influence of the wall and the quantum trajectory is clearly visible. 



In order to render the problem tractable with our 
present means, here we content ourselves with the as- 
sumption that a second absorbing wall is located at . 
This is a good approximation for a two-throat model in 
which the second throat is much "shallower" than the pri- 
mary inflationary throat (which is responsible for most 
of the exponential expansion). Our hope is that this toy 
model with absorbing walls at both 4)q and (/j^^ will give 
us an idea about the behavior to be expected in a more 
realistic case. Of course, placing a reflecting wall at (f)^^ 
is also technically possible but seems physically less jus- 
tified at present. Therefore, in the following, we concen- 
trate on the case of two absorbing walls. 

A. Geometrically Amended PDF 

Our goal is now to derive the PDF in presence of two 
absorbing barriers, i.e. the equivalent of Eq. (|87p but now 
with the two boundary conditions 

P(¥>-0o)--P(¥'-'^uv)=O. (88) 

There are two ways to derive the corresponding PDF. 
The first one is based on the method of images, see 



e.g. Refs. [H, [sj]. We start with the unrestricted 
PDF Pg(((f — 0mean), defined in Eq. ([551) . centered at 

mean 

= (pci + {5f2)- In order to "kiU" the contribution 
of Pg at (/>uv: introduce another "source" located at 

= 2(j)^^ - (/)mcan (this SOUrCC is the "image" of 0mcan 

with respect to the wall located at ) and weighed with 
a minus sign. In the same manner, in order to remove 
the contribution of Pg at 0o, we introduce a second image 

located at 0g^'' = 200 ~ ^moan still weighed with a minus 
sign (this source is the image of 0moan with respect to 
the wall located at 4>o). Therefore, the new PDF is given 
by 

Pg(^ - 0n.ca„) - Pg - ^'^^l) -P,{^- 0^'^) 
= Pg(</5 - 0,„oan) - Pg(<y9 - 2<j)^^ + (/)inoan) 

-Pg(</J-20o + 0mea„). (89) 

However, this new distribution function does not yet sat- 
isfy the boundary conditions ([55)) because the source 0^ 
now gives a contribution at (j)Q. In the same manner, the 
source (j)^^^ gives a contribution at (j)^^ . The cure is ob- 
viously to add extra images. Therefore, we introduce the 
source cf)^^^ which is the image of 4>l^^ with respect to the 
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wall located at 



I.e. 



^(2) _ fli ^ 

' TIV T UV 



(90) 



this process must be repeated ad infinitum^ i. e. one must 
introduce an infinite number of images. The location of 
nth source can be expressed as 



or 



200- 



(91) 



This time the PDF must be weighted with a plus sign 
since the new contribution is negative. We also introduce 
the image ^q'^'' , the image of (j)^^ with respect to the wall 
located at ^q. As a consequence, one obtains 



200. 



(92) 



But, as before, the new images will give new contribu- 
tions on the two walls. Clearly, in order to obtain a 
distribution that satisfies the boundary conditions 



0(2n) 
,(2n+l) 

i(2n+l) 
uv 



0moan + 2n (0o " 0uv ) i 
0mcan + 2n (0^^ ~ 4>o) , 

200 - 0mcan + 2n (0o " 0uv ) : 



20. 



2n\ 



(93) 
(94) 

(95) 
(96) 



Note the distinction between the positions of the even- 
and odd-numbered sources (n — 0,1,2...). Then, the 
complete distribution is obtained by taking the sum 
over all n contributions with their appropriate respec- 
tive signs, i.e. a minus for the "odd sources" and a plus 
for the "even sources" . The final expression is given by 



oo oo 



E 



exp ■ 



- 0n 



2n (0UV - 0o)] 



2 (^V'?) 



exp ■ 



[tf + 0nican - 200 + 2n (0„y - 0o)]^ 

2 (J^f ) 



(97) 



This formula is very general: it gives the PDF for any Gaussian process (characterized by the variance ((5</'i) ^^^id 
mean value {S1P2)) that is restricted between two absorbing barriers. 

As mentioned above, there exists another way to reach the above result. This second method consists in starting 
from the first line of Eq. (p8|) . i.e. the definition P = {5{(p — ip[£,])), but with a Dirac (5-function compatible with the 
boundary conditions of Eq. (j88p . The corresponding representation of the Dirac (5-function is given by 



■^sh 

71=1 



sm mr 



V' - 00 



(98) 



This is a J-function in the sense that /f"^ dipS{ip — ip') f{(p) — f{^') for any function f{(p) and any ip' e [0o, ' 



From this expression, it is straightforward to establish that 



00 

^exp 

n=l 


[ 2 (0 


00 

^exp 

n=l 


[ 2 (0 



nn 



-{(p-(Pc\ - {5(P2)) 



-{(p + (pcl + {S(p2) - 20o) 



(99) 



Then, using the identity 
g e-(^+"^)^ = ^-ff;^e-(^)%osf^), 

n— — 00 n— 1 ^ ' 

(100) 

it is easy to prove that Eq. is in fact exactly Eq. (1^71) . 
This PDF is the main result of this paper: it is appli- 
cable to any model where the range of variation of the 



stochastic inflaton field is limited by two absorbing walls. 
Finally, let us notice that the distribution (|97l) is not nor- 
malized as is expected since the two boundary conditions 
are absorbing walls. 
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B. Inflaton Between Two Absorbing Walls while the normahzation 3^ is given by 



In order to study the behavior of the inflaton field in 
a KS throat with its boundaries marked by absorbing 
walls, we now calculate the mean field value {tp) using 
the geometrically consistent PDF (|97)) obtained in the 
previous Section. It is given by the following expression 



1 



(101) 



where the "normalization" can be calculated from 



(102) 



00 



Lengthy but straightforward calculations show that 

n^n^ (Sifl) 



/•<Puv ^ 2 

/ P{ip)(pdip = > — exp 







2 {(t>vv - 4>o) 
cos(n7r) - (f)o] , (103) 
I 



J< = 



OO 

^ n 

n=l 

X sin 



■ exp 



niT (0n 



^ ■ (104) 



Finally, one can push this reasoning further to calculate 
the the variance of the field from (</'^) defined as 



1 rvvv 



dip P{p) p^ 



(105) 



00 



Using the same techniques, one obtains 



1 ^ 2 

n— 1 



2 (0UV - M'' 



X sm 



(106) 



— cos(ri7r) 



- 2 



It is interesting to check the consistency of the above 
expressions at the initial time (when ip — (/Sin)- Ini- 
tially, {5p2) = = the series in Eqs. (|103p 
and (|104p can be calculated explicitly using formu- 
las (1.441.1) and (1.441.3) of Ref. The resuh reads 
//o"^ d(/3P(iy9) 1^ = (^in and 3\f = 1. As a consequence, one 
checks that (</3)j,j = Pi-a as expected. 

The quantity {p>) is represented in Fig. [3] in the case 
where the two walls are respectively located at = 
lO^^TOpi and = lO^'^mpj. The infiuence of the walls 
is clearly visible from a comparison of the two curves. 
While the unbounded distribution artificially predicts 
that the brane can be outside the throat, the "true" quan- 
tum trajectory accounting for the geometric restrictions 
shows an interesting behavior: as the brane is approach- 
ing the wall, its position starts oscillating and then be- 
comes stabilized at a value > 4>o. Hence, the presence of 
the wall prevents if from violating its geometric limit <^o > 
i.e. from going beyond the bottom of the throat. 

Several words of caution are in order here. Firstly, 
the above conclusion is valid for particular initial con- 
ditions. Clearly, if the parameters are such that the 
stochastic correction {5ip2) remains small, the presence 



of the wall is not felt. Secondly, a much more serious 
problem is that the perturbative approach is not valid in 
the regime where the brane position is oscillating, i.e. on 
the left part of Fig. |31 Recall, however, that here we 
are considering the UV scenario in which the infiaton 
field value decreases during inflation. Therefore, Fig. [3] 
should be read "from right to left" : in the classical and 
the perturbative stochastic approach without the walls 
described by Eq. ([55]) (green dotted and red dashed lines, 
respectively), the brane travels towards the bottom of 
the throat, well within the regime of validity of our per- 
turbative treatment. Initially, also the trajectory found 
in the presence of the walls from Eq. ((W)) (black solid 
line) overlaps with these, but in the vicinity of 00, the 
"geometry-conscious" brane changes direction: it starts 
climbing upwards in the throat again, as can be seen in 
Fig. [31 This turnaround occurs within the (hatched) re- 
gion where the perturbative approach is reliable. 

In addition, it seems reasonable to conjecture that the 
oscillatory behavior mentioned above is real despite the 
fact that it occurs outside the regime of validity of our ap- 
proximation. The reason is as follows. In fact, the regime 
of validity indicates where the calculation of {5<p>\) and 
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{5lp2) is no longer reliable. But it does not limit in any 
way the validity of Eqs. and (IT(H| . Clearly, 

the oscillatory behavior comes from the peculiar struc- 
ture of the PDF in Eq. (P7|) which gives rise to the appear- 
ance of trigonometric functions in Eqs. (|103p and (|104p . 
As long as the stochastic effects grow (and no matter how 
quantitatively the do, i.e. no matter the detailed behavior 
of {Sipi) and {S(p2)), the brane will feel the wall at some 
point and, consequently, the trigonometric functions in 
Eqs. (|103p and (|104p will start to play a role. Hence, it 
is very likely that the brane position will oscillate even 
if, with the perturbative method used here, we cannot 
calculate the fine structure of these oscillations. 

Therefore, while the results obtained above probably 
can not be considered a fully realistic calculation of the 
quantum trajectory, we may well take them as an indi- 
cation that geometric limits are of great importance in 
stochastic DBI inflation. The finite size of the extra di- 
mensions, modelled by the presence of the two absorbing 
walls, changes the stochastic corrections to the classical 
field trajectory considerably. 

Finally, we have proven that, for the CKS potential 
with a constant term, stochastic effects can be dominant 
and occur near the bottom of the throat. The question of 
eternal inflation is clearly a more complicated issue. In 
particular, computing the stationary distribution in this 
case cannot be done using the present formalism, even if 
one can argue (as we have above) that the existence of 
an eternally inflating regime is likely in brane inflation. 



V. CONCLUSIONS 

We now conclude our investigation by revisiting our 
main results. In this paper, we have generalized the ap- 
proach of Refs. [IOjEI to DBI inflation models, in which 
the kinetic term of the inflaton is modified through an 
geometry-imposed upper limit on the field velocity. To 
this end, we solved the DBI Langevin equation using a 
perturbative expansion in the noise. It turns out that the 
results of Refs. [1^ [2l| essentially only change by addi- 
tional factors of 7, reflecting the fact that the distinction 
between long and short wavelength fluctuations, which 
is at the heart of the stochastic inflationary approach, 
now is deflned with respect to the sound (instead of the 
Hubble) radius scale. 

Our calculation yields easy expressions for the PDFs 
describing the probability for the patch-averaged inflaton 
to assume a given value in one Hubble domain, or in the 
entire Universe taking into account the size of each aver- 
aging domain. We calculated and plotted these PDFs for 
the example of chaotic Klebanov Strassler inflation, both 
for the cases with and without a constant term in the po- 
tential. In the absence of a constant term Vq, stochastic 
effects are found not to alter the classical field trajectory 
in a significant way. However, in both the UV (e = +1) 
and IR (e = —1) formulation of a potential including a 
constant term, the quantum behavior of the field plays a 



dominant role for small field values, i.e. at the bottom of 
the Klebanov Strassler throat. Hence, in the UV model, 
where the brane moves downwards in the throat, these 
effects occur at the end of infiation, whereas in the IR 
case, with the brane climbing up the throat, they come 
into play at the very beginning. 

The main result of this paper is our demonstration 
that an additional subtlety arises due to the geometric 
restriction which confines the field value of the coarse- 
grained field between (jjQ < tp < (j)^^ : to prevent the 
PDFs from extending below (f>Q (where the throat -and, 
consequently, the stringy spacetime itself- ends), we in- 
troduced two absorbing walls located at ipo and ^uv- 
This changes the PDF's shape, and within the bounds 
of validity it is the "absorbing walls probability density 
functions" that must be used to describe the stochastic 
behavior of the inflaton fleld. It should be interesting 
to put this technique to use for the "DBI version" of 
quintessence proposed in Ref. [46j . 

The question of existence of an eternally self- 
reproducing regime in the DBI case is more complicated. 
We have argued that such a regime does not exist if 
the potential is of the pure to^0^ type, but that in the 
case of an additional constant term Vq , eternal inflation 
seems possible. Given the intrinsic validity limitations 
of our present approach, a fully numerical treatment of 
the DBI Langevin equation would be necessary to answer 
this question. 

Finally, we would like comment on possible brane tra- 
jectories beyond the single throat model. Globally, the 
string geometric 6d background is given by a compacti- 
fled bulk Calabi Yau space (where the metric is unknown) 
whose "corners" may comprise several (possibly general- 
ized) Klebanov Strassler throats with different parame- 
ters. Some of these throats can contain additional anti- 
branes at their bottom, like in the original scenario of 
Ref. . One could therefore imagine that a brane starts 
out at the bottom of an IR throat, moving upwards to- 
wards the bulk, then transverses a short distance within 
the bulk before dropping down into a UV throat. 

Let us consider qualitatively the quantum effects ex- 
perienced by a brane following such a trajectory: While 
in the IR throat, the impact of stochastic inflation is 
strongest at the very bottom, i.e. quantum effects will 
push the brane upwards even faster than its classical evo- 
lution demands. Having made its way across the bulk, 
the brane will descend into the UV throat first on a purely 
classical trajectory. However, as it nears the bottom of 
the second throat, the stochastic infiuence grows again, 
pushing the brane away from the bottom and possibly 
upwards into the bulk again. At the limit, one may there- 
fore imagine a scenario where the brane is propelled out 
of a UV throat back into the bulk every time it comes 
sufficiently close to a throat's bottom. 

Since we ignore the formulation of the metric in the 
bulk, it seems difficult to carry out the full calculation 
for the brane trajectory sketched above. However, we do 
know how to describe the metric within each individual 
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throat. As a first toy model, we therefore considered the 
case of a brane trapped between two absorbing walls lo- 
cated at one throat's bottom and edge, respectively. This 
successfully prevents the brane from meandering "outside 
of spacetime" (below the bottom of the throat at (j)o), 
and we have therefore solved the severe conceptual issue 
stated in the introduction of this paper. Near the edge 
of the throat at (^^^ i the technical difficulty of describ- 
ing the stochastic effects remains because to-date the 6d 
bulk metric remains unknown. The next step would be to 
write down a joint warp factor for two throats of different 
depth which are "glued together" at their edges. Given 
the quantum behavior found in this paper, our expecta- 
tions might point us towards an increased probability to 
find the brane near the "throat matching point" . 



Lorentz factor 7 defined in Eq. ([6]) , we now introduce the 
additional hierarchy of DBI parameters Sj, the so-called 
"sound flow parameters", see e.g. Refs. [3, These 
are defined as 



Sn+l 



d In \S„ 
dN 



60 = 



7 



(A4) 



To find the "slow-roll" version of Eqs. Q and ([5]) (where 
by "slow-roll" we now understand that both of the pa- 
rameter sets ei,Si ^ 1), we use the DBI background 
equations and the definition of 7 in Eq. ([6]) to show that 
the following expression holds 



P ^—H^ei. 

K7 



(A5) 
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Appendix A: Slow-Roll Limit 

In this short appendix, we justify the use of the equa- 
tion ~ kV/3 in the DBI case, and present an alter- 
native derivation of the (classical) slow-roll DBI Klein 
Gordon equation. In the standard (non-DBI) case, a 
convenient tool to study the inflationary evolution is the 
following hierarchy of slow-roll parameters (47l - l49| 



dln|e^ 
dN 



If 



(Al) 



Then, straightforward manipulations of the background 
Einstein equations lead to the two following equations 



H' = 



K V 



3 - ei + -62 ) 



3 1 -ei/3 ' 



(A2) 



(A3) 



where a prime denotes a derivative with respect to the 
inflaton field. In the slow- roll limit where £1,62 <C 1, we 
therefore obtain the (Friedmann) equation ~ kV/3 
and the slow-roll version of the Klein Gordon equation, 
namely 3H4> ~ — F'. 

We now derive the same equations, but for DBI in- 
flation. In this case, the hierarchy (jAl[) is no longer 
sufficient and, in order to describe the evolution of the 



from which we can deduce that 



0(e2-2ei-(5i) 



(A6) 



To replace the terms proportional to T' in Eq. ([S]), we 
multiply the expression for Si (see Ref. [l3|) by and 
use Eq. (jA5|) to find 



T' 



V 

"7-1 



iH- 



7 



7-1 



H(t) 



7 



(72-l)(7-l) 



(A7) 



When we use these e^, 5i to rewrite Eqs. (|4]) and (|5|), we 
therefore find 



V 



£2 



3 1 
2ei 



2ei7/ (37 + 3) ' 

'^i , 3(7+1) 
(7 + 1) 



7 



7 + 1 V^' 
272 H 



(AS) 



(A9) 



As announced, in the limit where both £^,(5; ^ 1, one 
recovers that iJ^ ~ i^v/i and Eq. Note that in 

Eq. (|A9|) both the Ci and 5i are supposed small to obtain 
Eq. (ITUl) . while it sufficient to have £1 <C 1 to re-derive 
the slow-roll version of the Friedmann equation. 



Appendix B: Perturbed Klein Gordon Equation 

In this appendix, we briefly sketch how to obtain the 
formulation Eq. (ITTl) for the DBI Klein Gordon equa- 
tion at the linearly perturbed level. For the case of a 
standard inflaton field, it is well known that the Klein 
Gordon equation for linear perturbations of the field 5(j)k 
and scalar perturbations of the metric $fc reads [soj 

5'>Pk + 'iH54)k+ (% + V'] Sq^k = 4(/.$fe-2V^'$fc . (Bl) 
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One can express V' and V" in terms of the slow-roll pa- 
rameters defined in the Appendix El 

r ^4i^^(3-. + |), (B2) 

-2el^'j + le,e,-le,e,y (B3) 

where the function z = a(j)lH. Moreover, we have (soj 

^k^H€^-J4>u-H<^u. (B4) 

Therefore, using Eqs. (jB2P - ()B4p we can rewrite Eq. (|Bip 
as 

+ ^ eie2 - ^-^Hk = — <i>fc (2 - 2ei + e^) . (B5) 
z z / a 

Note that precisely the same equation is found by 
working backwards from the well-known evolution of 
the standard Mukhanov Sasaki variable [in confor- 
mal time ?7, where At = adr], Eq. (|B6p simplifies to 



v'^ + (fc2 - z"/z) Vk = 0, compare Eq. ((20)) ] . 

Vk + Hvk+i^ jvk = 0, (B6) 

which is defined as 

Vk ^ aS4>k + z^k, (B7) 

To see this, express the time derivatives of z in terms 
of the €i parameters, and use Eq. (IB4p as well as its 
derivative along with 

z k"^ ei-ff^ 

— ^*fc = —z— ie2 - 2ei) S<pk - ^iH S<j)k 
a a'' 2 

+eiH^-<^k,, (B8) 
a 

which is a consequence of the Einstein equations [50| . 

We now proceed to deriving the analogous equations 
for a scalar field with DBI dynamics. Not surprisingly, 
the expressions are significantly more involved in this 
case, and the perturbed DBI Klein Gordon equation 
reads 



r 



3 



6H ■ 3V' 3T' , 



7 



7 



7 



T 



— (37^ - 37^* - 47^ + 37^ + 37-2) 



a^7^ 



V" T", 3 ^ . 72 -IT' 



7' 
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27 T 



^^^6^^-^(37^-37^ + 2) 



(B9) 



Despite its complicated appearance, it is easy to see that 
this equation reduces to Eq. (jBl[) upon setting 7 = 1. 
Lengthy but straightforward calculations allow to express 
each term in Eq. (|B9[) in terms of the (6^,(5^) parameter 
hierarchy defined in Appendix [Xj in close analogy to the 
expressions (jB2l) and (jB3p in the standard case. 
In addition, the DBI version of Eq. (|B4[) reads 



ay 



3/2 



fc ■ 



(BIO) 



now with 



aj 



3/2, 



/H. Combining all of these ingre- 



dients, one arrives at Eq. ([TTl) in Section IIIB| which is 
equivalent to Eq. (IBSP for the case of standard inflation 
in which 7 = 1, (5^ = 0. 

Again, Eq. (|17p could also have been obtained from 



Vk + Hvk 



z + Hz 



(BU) 



where in the DBI case, too, the definition of the 
Mukhanov- Sasaki variable is 



Vk = a-i^l"^ 5(t)k + z^k 



(B12) 



but with the new z defined below Eq. (jBlOp . Here we 
have used 



z fc2 eiH'^ 

$fc = — — (£2 - 2ei - (5i)(5(/)fc - eii?(5(/>fc 



07^/2 0^7^ 



ay 



3/2 



(B13) 



as found from the Einstein equations with DBI scalar 
matter. 



20 



[1] L. McAllister and E. Silvcrstcin, Gen. Rel. Grav. 40, 565 [27 
(2008), ai-Xiv:0710.2951 [liep-th]. 

W. H. Kinney (2009), 0902.1529. [28 
D. Baumann and L. McAllister, Ann. Rev. Nucl. Part. 

Sci. 59, 67 (2009), 0901.0265. [29 
L. C. Lorcnz (2010), 1002.2087. [30 
G. R. Dvali and S. H. H. Tye, Phys. Lett. B450, 72 [31 
(1999), hep-ph/9812483. 

C. P. Burgess et al., JHEP 07, 047 (2001), hep- [32 
th/0105204. 

S. Kachru et al., JCAP 0310, 013 (2003), hep- [33 
th/0308055. 

D. Baumann ct al., JHEP 11, 031 (2006), hep- [34 
th/0607050. 

D. Baumann, A. Dymarsky, I. R. Klebanov, L. McAllis- 
ter, and P. J. Steinhardt, Phys. Rev. Lett. 99, 141601 [35 

(2007) , 0705.3837. 

D. Baumann, A. Dymarsky, I. R. Klebanov, and L. McAl- [36 
hster, JCAP 0801, 024 (2008), arXiv:0706.0360 [hep-th]. 
A. Krause and E. Pajer, JCAP 0807, 023 (2008), [37 
0705.4682. ' [38^ 

E. Pajer, JCAP 0804, 031 (2008), 0802.2916. 

D. Baumann, A. Dymarsky, S. Kachru, I. R. Klebanov, [39 
and L. McAllister, JHEP 03, 093 (2009), 0808.2811. 

E. Silvcrstein and D. Tong, Phys. Rev. D70, 103505 [40 
(2004), hep-th/0310221. 

M. Alishahiha, E. Silverstein, and D. Tong, Phys. Rev. [41 
D70, 123505 (2004), hep-th/0404084. 
L. Lorenz, J. Martin, and C. Ringeval, Phys. Rev. D78, 
063543 (2008), 0807.2414. 

L. Lorenz, J. Martin, and C. Ringeval, Phys. Rev. D78, [42 

083513 (2008), 0807.3037. 

A. A. Starobinsky (1986), in *De Vega, H.j. ( Ed.), 
Sanchez, N. ( Ed.): Field Theory, Quantum Gravity and [43 
Strings*, 107-126. 

A. A. Starobinsky and J. Yokoyama, Phys. Rev. D50, [44 
6357 (1994), astro-ph/9407016. [45 
J. Martin and M. Musso, Phys. Rev. D73, 043516 (2006), 
hep-th/0511214. [46 
J. Martin and M. Musso, Phys. Rev. D73, 043517 (2006), 
hep-th/0511292. [47 

F. Kuhnel and D. J. Schwarz, Phys. Rev. D78, 103501 

(2008) , 0805.1998. [48 
F. Finelh, G. Marozzi, A. A. Starobinsky, G. P. 
Vacca, and G. Venturi, Phys. Rev. D79, 044007 (2009), [49 
0808.1786. 

F. Kuhnel and D. J. Schwarz, Phys. Rev. D79, 044009 [50 

(2009) , 0810.5686. 

D. Seery, JCAP 0905, 021 (2009), 0903.2788. 
F. Kuhnel and D. J. Schwarz (2010), 1003.3014. 



F. Finelli, G. Marozzi, A. A. Starobinsky, G. P. Vacca, 
and G. Venturi (2010), 1003.1327. 

X. Chen, S. Sarangi, S. H. Henry Tye, and J. Xu, JCAP 
0611, 015 (2006), hep-th/0608082. 

A. Vilcnkin, Phys. Rev. D27, 2848 (1983). 

A. D. Linde, Phys. Lett. B175, 395 (1986). 

F. Helmer and S. Wiuitzki, Phys. Rev. D74, 063528 

(2006), gr-qc/0608019. 

A. J. Tolley and M. Wyman, JCAP 0804, 028 (2008), 
0801.1854. 

A. Sommerfeld, Electrodynamics: Lectures on Theoretical 

Physics (Academic Press, New York, US, 1964). 
A. Sommerfeld, Partial Differential Equations in 
Physics: Lectures on Theoretical Physics (Academic 
Press, New York, US, 1964). 

J. Garriga and V. F. Mukhanov, Phys. Lett. B458, 219 

(1999), hep-th/9904176. 

X. Chen, Phys. Rev. D71, 063506 (2005), hep- 
th/0408084. 

X. Chen, JHEP 08, 045 (2005), hep-th/0501184. 

R. Bean, J. Dunkley, and E. Pierpaoli, Phys. Rev. D74, 

063503 (2006), astro-ph/0606685. 

R. Bean, X. Chen, H. V. Peiris, and J. Xu, Phys. Rev. 
D77, 023527 (2008), 0710.1812. 

L. Lorenz, J. Martin, and C. Ringeval, JCAP 0804, 001 
(2008), 0709.3758. 

M. Abramowitz and L A. Stegun, Handbook of mathe- 
matical functions with formulas, graphs, and mathemat- 
ical tables (National Bureau of Standards, Washington, 
US, 1970), ninth ed. 

I. S. Gradshteyn and I. M. Ryzliik, Table of Integrals, 
Series, and Products (Academic Press, New York and 
London, 1965). 

S. Gratton and N. Turok, Phys. Rev. D72, 043507 
(2005), hep-th/0503063. 

S. Chandrasekhar, Rev. Mod. Phys. 15, 1 (1943). 

J. Martin and M. A. Musso, Phys. Rev. D71, 063514 

(2005), astro-ph/0410190. 

J. Martin and M. Yamaguchi, Phys. Rev. D77, 123508 
(2008), 0801.3375. 

D. J. Schwarz, C. A. Terrero-Escalante, and A. A. Garcia, 

Phys. Lett. B517, 243 (2001), astro-ph/0106020. 

S. M. Leach, A. R. Liddle, J. Martin, and D. J. Schwarz, 

Phys. Rev. D66, 023515 (2002), astro-ph/0202094. 

D. J. Schwarz and C. A. Terrero-Escalante, JCAP 0408, 

003 (2004), hep-ph/0403129. 

V. F. Mukhanov, H. A. Feldman, and R. H. Branden- 
berger, Phys. Rept. 215, 203 (1992). 



